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ABSTRACT 


New  labeling  techniques  are  provided  for  accelerating  the  basis 
exchange  step  of  specialized  linear  programming  methods  for  network 
problems.  Computational  results  are  presented  which  show  that  these 
techniques  substantially  reduce  the  amount  of  computation  involved  in 


updating  operations. 


1.0  INTRODUCTION 


In  solving  minimum  cost  flow  network  problems  by  specialized  linear 
programming  methods,  an  important  question  is:  How  can  one  update  the 

spanning  tree  basis  with  the  least  amount  of  effort?  A partial  answer  to 
this  question  is  provided  by  special  list  structure  techniques  such  as 
the  API  method  [6]  and  the  more  recent  ATI  method  [9],  which  have  con- 
tributed dramatically  to  improving  the  efficiency  of  network  algorithms 
(see,  e.g.,  [3,  5,  6,  7,  11]).  This  paper  addresses  the  issue  of  which 
supplemental  techniques  can  be  used  to  enable  these  list  structures  (and 
particularly  the  ATI  method)  to  be  implemented  with  greater  efficiency. 

As  shown  in  [6],  the  major  updating  calculations  of  a basis  exchange 
step  can  be  restricted  to  just  one  of  the  two  subtrees  created  by  dropping 
the  outgoing  arc.  Consequently,  a natural  goal  is  to  identify  the  smaller 
of  these  two  subtrees  by  means  of  a function  t(x)  that  names  the  number 
of  nodes  in  the  subtree  "headed  by  node  x."  A clever  and  rather  intricate 
procedure  for  doing  this  was  proposed  by  Srinivasan  and  Thompson  in  ] 1 3 1 . 
Unfortunately,  this  procedure  requires  sorting  the  nodes  of  the  subtree 
by  their  distances  from  the  root,  and  then  further  entails  a full  subtree 
update  of  both  the  distance  values  and  the  t(x)  values  at  each  basis 
exchange  step.  Because  of  the  substantial  amount  of  work  required  to 
update  the  t(x)  list,  the  advantages  of  using  this  list  have  been  largely 
offset  by  the  computational  costs  involved  in  its  maintenance,  and  the 
potential  of  the  original  Srinivasan-Thompson  proposal  has  not  been  fully 
realized. 
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The  purpose  of  this  paper  is  to  propose  a new  type  of  relabeling 
scheme  that  succeeds  in  updating  t(x)  without  sorting.  In  fact,  this 
scheme  requires  even  less  work  than  to  update  the  distance  values  of  [13]. 
The  relabeling  is  based  on  "absorbing"  t(x)  into  the  updating  calcula- 
tions of  the  ATI  method.  Moreover,  these  calculations  are  carried  out 
simultaneously  with  the  procedures  [9]  for  updating  other  changes  intro- 
duced by  the  basis  exchange  step. 

To  achieve  the  integration  of  the  ATI  calculations  and  the  update  of 
t(x),  an  index  function  f(x)  is  introduced  that  names  the  last  node  in 
the  subtree  rooted  at  x.  Additionally,  it  is  shown  that  f(x)  makes  it 
possible  to  streamline  the  ATI  calculations.  Finally,  as  a bonus,  it  is 
shown  that  t(x)  can  accommodate  all  of  the  relevant  functions  filled  by 
the  distance  values,  and  hence  can  replace  these  values.  Computational 
results  presented  in  the  last  section  indicate  that  the  net  gains  of  all 
these  advantages  produce  a substantially  improved  procedure  for  imple- 
menting the  basis  exchange  operations. 

2.0  NOTATION  AND  DEFINITIONS 

It  will  be  assumed  that  a basis  tree  with  n nodes  and  n - 1 arcs  is 
known  and  has  been  rooted.  The  root  node  will  be  regarded  as  the  "highest" 
node  in  the  rooted  tree  with  all  other  nodes  hanging  below  it.  If  nodes 
i and  j denote  endpoints  of  a common  arc  in  the  rooted  tree  such  that 
node  i is  closest  to  the  root,  then  i is  called  the  predecessor  of  node  j 
and  node  j is  called  an  immediate  successor  of  node  i. 

The  following  notational  conventions  will  be  used  to  identify  the 
components  of  the  basis  exchange  step: 
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(p,q)  = the  arc  leaving  the  basis,  where  p is  currently  the  predecessor 

of  q. 

(u,v)  = the  arc  entering  the  basis,  where  u is  the  node  whose  unique 
path  to  the  root  node  contains  arc  (p,q). 

T = the  basis  tree. 

T(x)  = the  subtree  of  T that  is  rooted  at  node  x (hence  the  subtree 
that  includes  x and  all  its  successors  under  the  predecessor  ordering). 

p(x)  = the  predecessor  of  node  x where  p(x)  = 0 if  node  x is  the 
root  node. 

. s(x)  = the  "thread  successor"  of  x. 

Intuitively,  function  s may  be  thought  of  as  a thread  which  passes 
through  each  node  exactly  once  in  a top  to  bottom,  left  to  right  order 
starting  from  the  root  node. 

More  precisely  the  function  s satisfies  the  following  inductive  charac- 
teristics : 

(a)  Letting  1 denote  the  root  node,  the  set  { l , s (1 ) , s2 ( 1 ) , . . . sP  ^ ( 1 ) } 

j 

is  precisely  the  set  of  nodes  of  the  rooted  tree  where  s2(l)  = s(s(l)), 
s'*  = s(s2(l)),  etc.  The  nodes  1,  s(l),...sk  *(1),  will  be  called  the 

k 

antecedents  of  node  s (1). 

(b)  For  each  node  i other  than  node  s°  ^(1),  s(i)  is  one  of  the 
immediate  successors  of  node  i,  if  i has  a successor.  Otherwise,  s(i)  is 

» 

an  immediate  successor  of  the  closest  predecessor  of  node  i,  say  node  x, 

i 

such  that  node  x has  an  immediate  successor  which  is  not  an  antecedent  of 
node  i. 

(c)  sn(l)  » 1;  that  is,  the  last  node  of  the  tree  threads  back  to  the 

i 

root  node. 
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By  virtue  of  the  foregoing  characterization,  the  set  of  nodes  of 
1 k 

T(x)  is  {x,s  (x),...,s  (x)}  where  k is  the  largest  number  such  that 

k ] k-1 

p(s  (x))  is  one  of  the  nodes  x,  s (x),...,s  (x) . By  convention 

sX(x)  = s(x)  and  s^(x)  = x. 

t(x)  = the  number  of  nodes  in  T(x). 

j* 

f(x)  = the  "last  node,"  s (x) , of  the  nodes  in  T(x),  where  r = t(x)  - 1. 

Figure  1 illustrates  the  above  functions  as  follows.  The  NODE  array 
specifies  the  node  names.  The  entries  in  the  arrays  p,  s,  f,  and  t parallel 
to  a node  name  specify  the  values  of  the  functions  p,  s,  f,  and  t for  that 
node  name.  Note  that  the  direction  of  the  arcs  in  Figure  1 corresponds  to 
the  direction  induced  by  the  predecessor  ordering  and  may  not  correspond  to 
the  direction  of  the  arc  in  the  underlying  problem. 

The  basis  exchange  step  may  be  visualized  as  consisting  of  two  components 

(1)  Dropping  arc  (p,q)  to  create  two  independent  subtrees:  T(q)  and 

T-T(q)  (where  the  latter  is  the  subtree  of  T that  excludes  T(q)  and  all  its 
nodes,  and  hence  which  excludes  the  "connecting  arc"  (p,q)); 

(2)  Adding  arc  (u,v)  to  create  a single  new  basis  tree. 

The  subtrees  T(q)  and  T-T(q)  can  be  viewed  as  any  two  node-disjoint 
trees  which  are  to  be  joined  by  an  arc  to  create  a new  basis  tree.  Updating 
operations  will  be  developed  first  to  make  T(q)  and  T-T(q)  into  label 
"independent"  trees  in  preparation  for  selecting  which  is  to  be  the  new 
"upper  tree"  (called  T^)  rooted  at  x^  and  which  is  to  be  the  new  "lower 
tree"  (called  T^)  rooted  at  It  may  then  be  assumed  that  the  root  of 

becomes  the  root  of  the  new  basis  tree.  Additionally,  (y^*y£)  will  be 
used  to  denote  the  arc  that  joins  T^  and  T2  where  y^  is  a node  of  T^  and 
is  a node  of  T2.  Next,  operations  will  be  developed  to  re-root  T2  at 
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in  preparation  for  attaching  7^  to  via  arc  (y^.y^)  to  create  the  new 
basis  tree.  (There  is  no  requirement  that  be  distinct  from  or  that 
y2  be  distinct  from  . ) 


3 . 0 UPDATING  OPERATIONS 

Using  the  preceding  definitions,  rules  will  be  given  to  find  updated 
values  p*(x),  s*(x),  t*(x),  f*(x)  for  p(x),  s(x),  t(x),  and  f(x). 

I . Update  p(x),  s(x),  t(x),  and  f(x)  to  make  the  subtrees  T(q)  and 
T-T(q)  independent. 

The  reader  may  find  it  helpful  to  perform  the  following  operations 
using  Figure  1. 

1.1.  Update  for  T-T(q) : 

For  p*(x) : No  updating  of  any  p(x)  is  required. 

For  s*(x) : Identify  the  node  y in  T-T(q)  such  that  s(y)  = q. 

Then  set  s*(y)  = s(f(q)).  No  other  s(x)  values  are  changed.  (Note:  While 

this  step  is  obviously  facilitated  by  directly  accessing  f(q),  the  identifi- 
cation of  y is  further  speeded  by  utilizing  the  f(x)  function  as  follows: 
First,  let  y1  = p.  Second,  if  s(y1)  = q,  then  y = y1  and  the  process  stops. 
Otherwise,  let  y^  = f(s(y^))  and  repeat  the  second  step.) 

For  t*(x):  Set  t*(x)  = t(x)  - t(q)  for  those  nodes  x on  the  path 

from  p to  the  root  of  T.  It  is  important  to  note  that  due  to  cancellation 
effects  of  subsequent  calculations,  this  step  can  be  restricted  to  the 
partial  backward  path  from  p to  the  "intersection"  node  z that  is  the  unique 
node  of  the  basis  loop  that  lies  on  the  predecessor  paths  from  both  u and  v, 
to  the  root  of  T,  excluding  node  z itself  from  consideration. 

For  f*(x) : Let  x*  = p(s(f(q))>.  If  p(s(f(q)))  = 0,  then  set  x* 


equal  to  the  root  node.  Set  f*(x)  = y (i.e.,  the  node  y determined  in 
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updating  s above)  for  those  nodes  x on  the  path  from  p to  x*,  excluding  x* 
itself  if  p(s(f(q)))  t 0.  (If  x*  = p,  then  no  updating  Ls  done  unless  x* 
is  the  root  node.) 

1.2.  Update  for  T(q) : 

For  p*(x) : Set  p(q)  = 0 

For  s*(x) : Set  s*(f(q))  = q. 

No  other  updating  of  any  p(x),  s(x),  t(x),  or  f(x)  values  is 
required  for  T(q) . 

I I . Decide  which  of  T(q)  and  T-T(q)  is  to  be  Tp  and  which  is  to  be  T2 . 

The  rule  to  use  for  picking  which  subtree  to  re-root  heavily  depends  on 
when  the  dual  variables  are  updated.  If  the  updating  of  the  dual  variables 
is  not  integrated  with  the  other  updating  operations,  it  will  be  clear  from 
the  operations  performed  in  Step  III  that  the  subtree  to  re-root  should  be 
selected  according  to  the  number  of  nodes  in  the  predecessor  paths  from 
v to  the  root  of  T-T(q)  and  from  u to  q.  In  particular,  the  subtree  asso- 
ciated with  the  path  containing  the  minimum  number  of  nodes  should  be  re- 
rooted. Further,  if  the  updating  operations  are  not  integrated,  the  dual 
variables  should  be  updated  in  the  subtree  containing  the  smallest  number 
of  nodes. 

Thus  the  following  procedures  will  be  computationally  investigated  in 
section  5: 

Procedure  1:  Let  T^  be  the  subtree  whose  predecessor  path  from  the 

proposed  new  root  to  the  old  root  is  smallest  and  separately  update  the  dual 
variables  in  the  subtree  containing  the  fewer  number  of  nodes. 

Procedure  2:  Same  as  procedure  1 except  that  if  T^  is  to  be  botli  re- 

rooted and  its  dual  variables  updated,  then  integrate  the  updating  of  the 
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dual  variables  with  the  other  updating  operations. 

One  computational  disadvantage  of  the  above  procedures  is  that  the 
number  of  nodes  in  the  path  from  u to  the  root  of  T-T(q)  is  not  known  and 
to  calculate  it  involves  traversing  the  path  from  the  intersection  node  v. 
to  the  root  of  T-T(q)  for  no  purpose  other  than  to  calculate  this  number. 
Another  computational  disadvantage  simply  involves  the  fact  that  there  are 
certain  computational  advantages  to  always  performing  all  updating  opera- 
tions on  one  subtree.  To  partially  overcome  these  difficulties,  a compro- 
mise procedure  is  proposed: 

Procedure  3:  If  t(q)  exceeds  n/2  let  T(q)  be  T and  let  T-T(q)  be 
(hence  x^  = q and  x^  = the  root  of  T).  Otherwise,  let  T-T(q)  be  and 
let  T(q)  be  T^  • Perform  all  updating  operations  on  integrating  the  dual 
variable  updating  with  the  other  updating  operations. 

Procedure  4:  Let  T^  be  the  subtree  containing  the  smallest  number  of 

nodes . 

III.  Make  )>2  the  new  root  of  1’2  and  update  p (x)  , s (x)  , t (x) , and  f (x ) . 

It  is  possible  to  combine  the  updating  of  the  functions  in  several 
computationally  efficient  ways,  for  it  is  possible  to  simultaneously  update 
all  of  the  functions,  while  making  exactly  one  traversal  of  the  nodes  in 
the  unique  predecessor  path  from  y^  to  x;  and  their  descendants.  However, 
for  readability  the  updating  of  each  function  will  be  presented  separately 
assuming  that  no  other  functions  have  been  updated  since  step  I.  Thus  t he* 
order  of  the  updating  procedures  given  below  is  arbitrary.  In  the  following. 
p(x),  f(x),  t(x),  and  s(x)  refer  to  the  value  of  the  functions  after  the 
execution  of  step  I.  (Again,  the  reader  may  find  it  helpful  to  carry  out 
these  operations  using  an  updated  Figure  1 and  procedure  3 from  section  II.) 
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For  p* (x) : Reverse  the  predecessor  orientation  of  the  path  from 

to  That  is,  if  x.(  = y0,  do  nothing.  Otherwise,  recursively  set 

P*(p(y))  = y and  then  y = p(y)  until  and  including  p(y)  = x„  initially 
setting  y = y Finally  set  p*(y.>)  = 0.  (It  is  important  to  note  that 
p(x)  and  p*(x)  must  be  kept  distinct  from  each  other;  i.e.,  it  is  not 
legitimate  to  replace  p(x)  by  p*(x)  in  the  computation.)  No  other  updating 
of  p(x)  is  required. 

For  t*  (x)  : Set  t*(y->)  = t(x^).  (But  if  T(q)  = , and  the 

restricted  update  of  t(x)  was  carried  out  in  step  1.1,  set  t*(y2>  = 
t(x  ) - t (q)  . ) Then  for  each  x on  the  predecessor  path  from  y„  to  x 
(excluding  x = x,,),  set  t*(p(x))  = t*(y2)  - t(x).  (Again  note  that  t(x) 
and  t*(x)  must  be  kept  distinct  and  it  is  assumed  that  the  function  p(x) 
has  not  been  updated.)  No  other  updating  of  t(x)  is  required. 

For  s*(x) : If  x^  = no  updating  of  s(x)  is  required.  Otherwise 

execute  the  following: 

i)  Set  x = y^,  w = f(x),  and  z = s(w). 

ii)  Find  the  unique  node  y such  that  s(y)  = x.  (Note  that  y may  be 
found  as  in  step  1.1  by  letting  y1  = p(x).) 

iii)  If  p(z)  = p (x ) , set  s*(y)  = z,  s*(w)  = p(x),  w = f(p(x)), 
z = s(w),  and  go  to  step  iv.  Otherwise  set  s*(w)  = p(x), 
w = y end  proceed, 

iv)  Set  x = p (x) . 

v)  If  x i4  x^  go  to  step  ii.  Otherwise  set  s*(w)  = y ^ and  stop. 

In  the  updating  of  f,  this  last  value  of  w (i.e.,  the  node 
whose  new  thread  is  the  new  root  y^)  plays  a primary  role. 

(Note  that  w is  either  equal  to  f(x^)  or  the  last  y.) 


For  f*(x):  If  x,  = y9  no  updating  is  required.  Otherwise,  let 

x^  be  the  unique  node  on  the  predecessor  path  from  to  x^  such  that 
p(x2)  = x^.  If  f(x^)  ^ f(x2)  then  f*(x2>  = f(x2)  (i.e.,  f(x2)  does  not 
change).  Otherwise,  f*(x^)  = y where  y is  the  unique  node  such  that 
s(y)  = x2-  (This  y may  be  found  as  in  I . 1 ; but  as  noted  above,  it  will  be 
identified  automatically  at  the  conclusion  of  updating  s.)  In  either  case, 
set  f*(x)  = f*(x2)  for  all  x on  the  predecessor  path  from  y2  to  x2-  No 
other  changes  are  made. 

It  appears  that  in  order  to  optimize  computational  operations,  step 
III  should  be  implemented  as  follows.  One  should  traverse  the  path  from 
y2  to  x2  simultaneously  updating  the  functions  p(x),  s(x),  and  t(x),  the 
basic  flow  values  associated  with  arcs  on  this  path,  and  co-ordinating 
the  updating  of  the  node  potential  values  with  the  updating  of  s(x). 

Next  the  new  predecessor  path  from  x2  to  y2  should  be  traversed  to  update 

f(x). 

IV.  Attach  to  by  adding  arc  (y],y2)  to  create  the  new  basis  tree 
(where  T2  is  now  rooted  at  y2  as  a result  of  step  III). 

As  before,  in  the  following  p(x),  s(x),  f(x),  and  t(x)  refer  to  the 
value  of  the  functions  p,  s,  f,  and  t after  the  execution  of  step  III  above 
and  the  rules  for  updating  each  function  assume  that  no  other  functions 
have  been  updated  since  step  III. 

For  p*(x) : Set  p*(y2)  = y^ 

For  s*(x) : Set  s*(f(y2))  = s(yj)  and  s*(yl)  = y2> 

For  t*(x) : Set  t*(x)  = t(x)  + t(y9)  for  all  x on  the  path  from 


y^  to  Xj.  (But  if  T(q)  = 1'2,  and  the  restricted  update  of  t(x)  was  applied 
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in  step  1.1,  then  the  current  step  should  be  restricted  to  those  x on  tin* 
path  from  y^  to  z,  excluding  z itself,  lor  the  "intersect  ion"  node  z as 
identified  in  1.1.) 

For  f * (x) : Set  x = p(s(y^)).  If  x = 0 set  x = Xj . Then  for 
those  nodes  x on  the  backward  path  from  y^  to  x,  excluding  x itself  if 
ptsCy^)  ^ 0,  set  f*(x)  = f(y2>. 

(Figure  2 indicates  the  results  of  applying  steps  1,  II,  III,  and  IV 
above  on  Figure  1.) 

The  proposed  procedure  for  updating  t(x)  clearly  requires  less 

effort  than  updating  the  distance  function  of  [13],  which  involves  an 

addition  for  every  node  of  the  subtree  T(q),  and  in  the  case  of  T(q)  = T , 

requires  an  addition  for  every  node  of  T.  The  fact  that  t(x)  can  replace 

the  distance  function  is  a direct  consequence  of  the  observation  that  t(x) 

can  be  used  in  essentially  the  same  manner  as  the  distance  function  to 

facilitate  operations  of  identifying  the  loop  created  by  adding  arc  (u,v) 

to  the  basis  tree.  Specifically  this  loop  can  be  located  as  follows: 

i)  Set  x = u and  x = v. 
u v 

ii)  If  t(x  ) < t(x  ) go  to  step  v. 
u v 

iii)  If  t(x  ) > t (x  ) go  to  step  vi . 
u v 

iv)  If  x^  t x^,  go  to  step  v.  Otherwise  stop.  The  loop  created 

by  adding  (u,v)  has  been  traversed  and  x = x = z is  the  unique  node  of 

u v 

this  loop  which  lies  on  the  predecessor  paths  from  both  u and  v to  the 
root  node  of  the  tree. 

v)  Search  the  predecessor  path  of  u starting  at  x^  for  a node  x 
such  that  t(x)  > t(xv>  and  set  xy  = x.  If  t(x^)  = t(xy)  go  to  step  iv. 
Otherwise  go  to  step  vi . 
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vi)  Search  the  predecessor  path  of  v starting  at  for  a node  x 

such  that  t(x)  > t(x  ) and  set  x = x.  If  t(x  ) = t(x  ) go  to  iv. 

u v u v 

Otherwise  go  to  step  v. 


4 . 0 INITIALIZATION 

It  is  left  to  characterize  the  procedure  for  establishing  the  initial 

values  of  t(x)  and  f(x).  This  occurs  simultaneously  with  the  initial 

determination  of  the  s(x)  values  as  follows. 

Let  x denote  the  root  of  the  tree.  Consider  the  step  in  which 
o ' 

k+1  k k 

s (x  ) = s(s  (x  ))  is  identified  (k  > 0).  If  s (x  ) is  the  predecessor 
o o o 

k+1 

of  s (xq)  (via  the  predecessor  indexing),  do  nothing.  Otherwise,  for  all 

i k 

nodes  s (x  ) on  the  predecessor  path  from  s (x  ) to  the  predecessor  of 
o o 

i 

s (x  ) , excluding  the  predecessor  of  s (x  ) itself,  set  t(s  (x  ))  = 
o o o 

k + 1 - i,  and  set  fCs^x  ))  = s^(x  ). 

o o 

When  the  last  node  s°  ^(x^)  of  the  network  is  determined,  set 

t(s^(x  ))  = n - i and  set  f(s*(x  ))  = sn  \x  ) for  all  s*(x  ) on  the  pre- 
o o o o 

decessor  path  from  s'1  * (x  ) to  x . 

o o 

To  easily  keep  track  of  the  index  i for  each  node  sX(x  ) that  is  to 

o 

be  considered  on  a given  step,  it  is  convenient  to  keep  a list  that  con- 


sists precisely  of  the  indexes  i of  the  nodes  s (x  ) to  x . Specifically, 

o o 

to  begin  with  the  list  contains  the  single  index  0 (for  s^(x  )).  When 

o 

k+1 

s (x  ) is  created,  the  number  k + 1 is  added  to  the  end  of  the  list, 
o 

When  a predecessor  path  from  s (xo>  is  traced,  consisting  of  r nodes  (say) 

s*(x  ) whose  values  t(s*(x  ))  and  f(s*(x  ))  are  to  be  set,  the  indexes  of 
o o o 


these  r nodes  will  be  exactly  the  corresponding  last  r numbers  on  the  list. 


13 


By  removing  these  numbers  from  the  list  just  before  adding  the  number 
k + 1,  the  desired  structure  of  the  list  is  maintained. 

5 . 0 COMPUTATIONAL  ANALYSIS 

5.1  Highlights  of  the  Development  of  the  ARC- 1 1 Computer  Code 

To  evaluate  the  foregoing  procedures,  henceforth  referred  to  as  the 
Extended  Threaded  Index  (XT1)  Method,  we  developed  a new  in-core  computer 
code  entitled  ARC-II  for  solving  capacitated  transshipment  problems. 

ARC-II  is  written  in  a "manilla"  FORTRAN  with  several  subroutines  to 
perform  the  various  updating  operations,  for  the  following  reasons: 

(1)  this  modular  approach  simplifies  testing  of  different  updating  pro- 
cedures, (2)  minimal  recoding  is  required  to  fit  different  machine  and 
computer  conventions,  and  (3)  unbiased  comparisons  can  be  made  with  codes 
which  have  not  been  "customized"  to  a particular  machine  or  compiler. 

One  disadvantage  of  this  approach,  of  course,  is  that  the  reported  times 
are  conservative,  since  programs  which  have  been  "tuned"  to  a particular 
operating  environment  execute  substantially  faster.  However,  our  purpose 
in  developing  ARC-II  was  to  obtain  unbiased  comparisons  between  the  XTI 
approach  and  other  procedures  available  in  the  literature.  To  this  end, 
the  same  starting  and  pivoting  procedures  as  described  in  [8,  9]  for  trans- 
shipment problems  are  used. 

After  initially  developing  ARC-II,  preliminary  testing  was  conducted 
on  the  recoding  rules  described  in  part  II  of  section  3.  Our  initial 
testing  indicated  that  procedure  3 was  never  a good  rule  and  that  procedure 
A marginally  dominated  procedures  1 and  2.  Thus,  the  times  on  the  ARC-II 
code  reported  subsequently  in  this  section  reflect  the  use  of  procedure  A. 
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Another  supplementary  feature  tested  was  an  "inverse  thread  function." 
Using  this  function,  in  conjunction  with  the  other  functions  previously 
discussed,  one  can  eliminate  all  searching  involved  in  the  basis  exchange 
operation;  i.e.,  updating  the  thread  function  is  simply  a matter  of 
resetting  known  pointers.  The  disadvantages  of  using  the  inverse  thread 
include  increased  memory  requirements  and  the  need  to  maintain  an  addi- 
tional set  of  function  values.  Our  tests  using  the  inverse  thread  function 
indicate  that  solution  times  are  reduced  by  approximately  5%.  In  our 
opinion,  this  reduction  is  not  sufficient  to  warrant  the  utilization  of 
additional  memory  space,  and  therefore,  ARC-II  does  not  use  such  a 
function. 

5.2  Computational  Comparisons 

To  determine  the  efficiency  of  the  XT1  procedures,  ARC-II  was  compared 
with  three  out-of-kilter  codes  referred  to  hereinafter  as  SHARE  [4],  SUPERK 
[1],  and  BSRL  (developed  by  T.  Bray  and  C.  Witzgall  while  at  Boeing 
Scientific  Laboratories).  Additionally,  one  dual  simplex  based  code  [5], 
called  DNET,  one  negative  cycle  code  [2],  called  BENN,  and  two  primal 
simplex  based  codes,  called  PNET,  PNET-I,  were  tested  for  comparative 
purposes.  PNET  [8}  uses  the  API  list  structure  [6]  and  PNET-I  [9]  uses 
the  ATI  list  structure  [9]. 

All  of  the  above  mentioned  codes  are  in-core  codes;  i.e.,  the  program 
and  all  of  the  problem  data  simultaneously  reside  in  fast-access  memory. 

All  are  coded  in  FORTRAN  and  none  of  them  (including  ARC-II)  have  been 
tuned  (optimized)  for  a particular  compiler.  All  of  the  problems  were 
solved  on  the  CDC  6600  at  the  University  of  Texas  Computation  Center  using 
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the  RUN  compiler.  The  computer  jobs  were  executed  during  periods  when 
the  machine  load  was  approximately  the  same,  and  all  solution  times  are 
exclusive  of  input  and  output;  i.e.,  the  total  time  spent  solving  the 
problem  was  recorded  by  calling  a Real  Time  Clock  upon  starting  to  solve 
the  problem  and  again  when  the  solution  was  obtained. 

Since  the  t 'St  problems  of  [11]  are  currently  used  worldwide  for 
comparison  purposes,  they  were  also  used  in  our  comparison.  This  com- 
parison included  several  different  types  of  problems  (e.g.,  assignment, 
transportation,  and  minimum  cost  flow  network  problems),  both  capacitated 
and  uncapacitated,  and  with  varying  node  and  arc  requirements.  The  pro- 
blem specifications  of  these  35  problems  as  required  on  the  input  cards 
to  the  network  generator  are  given  in  [11].  Problems  1-5  are  100  x 100 
transportation  problems;  problems  6-10  are  150  x 150  transportation  pro- 
blems. Problems  11-15  are  200  x 200  assignment  problems.  Problems  16-27 
are  400  node  capacitated  transshipment  problems;  problems  28-35  are  un- 
capacitated 1000  and  1500  node  transshipment  problems.  Table  1 reports 
the  solution  times  for  each  of  these  problems  for  each  of  the  codes. 

The  results  in  Table  I clearly  indicate  the  superiority  of  ARC-II 
over  all  other  codes  tested.  Furthermore,  the  data  indicate,  rather 
startlingly,  that  ARC-11  is  approximately  twice  as  fast  as  one  of  the  (pre- 
viously) fastest  codes  in  the  literature,  PNET-I. 

It  is  also  particularly  noteworthy  that  the  solution  times  for  the 
ARC-11  code  are  based  on  using  the  simple  pivot  strategies  of  [8],  rather 
than  the  more  sophisticated  candidate  list  strategies  whose  superiority 
has  been  documented  by  Mulvey  [12].  We  have  used  the  simpler  pivot 
strategies  to  make  it  possible  to  more  clearly  differentiate  the  contribution 
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TABLE  I 

Solution  Times  in  Seconds  on  a CDC  6600 


Problems 

ARC- 11 

PNET 

PNET-I 

DNET 

BENN 

SUPERK 

SHARE 

BSRL 

1 

.78 

1.30 

1.07 

12.85 

20.25 

5.68 

17.76 

30.25 

2 

.89 

1.49 

1.25 

13.56 

24.36 

6.47 

21.34 

21.59 

3 

1.01 

1.94 

1.64 

21.44 

34.56 

6.87 

26.16 

31  .47 

4 

.95 

1.64 

1.27 

17.96 

31.45 

6.57 

25.13 

36.47 

5 

1.25 

1.88 

1.63 

23.34 

52.10 

6.77 

30.97 

47.73 

6 

2.11 

3.55 

2.86 

46.10 

61.00 

11.05 

46.40 

46.64 

7 

2.23 

4.06 

3.37 

74.88 

DNR 

12.86 

65.92 

113.12 

8 

2.99 

4.72 

4.10 

97.92 

DNR 

13.69 

81.00 

175.10 

9 

2.99 

4.80 

4.15 

101.65 

DNR 

13.40 

81.21 

186.99 

10 

4.02 

5.88 

5.27 

95.96 

DNR 

14.13 

84.24 

184.75 

11 

1.92 

3.52 

2.31 

19.87 

17.44 

6.44 

19.93 

30.39 

12 

2.36 

4.87 

3.71 

26.58 

20.31 

6.47 

21.17 

22.08 

13 

3.13 

5.52 

3.47 

27.98 

24.92 

7.25 

25.81 

20.02 

14 

2.96 

6.02 

3.44 

30.15 

27.40 

6.95 

24.95 

23.11 

15 

3.12 

6.50 

4.79 

31.57 

DNR 

7.56 

27.05 

21.08 

16 

1.38 

_ . 40 

2.15 

14.77 

11.77 

5.27 

21.51 

15.05 

17 

1.87 

3.11 

2.60 

DNR 

20.10 

8.36 

32.40 

64.64 

18 

1.26 

1.92 

1.70 

DNR 

11.31 

5.13 

20.06 

18.31 

19 

1.72 

2.60 

2.40 

DNR 

20.62 

8.49 

31.75 

61.07 

20 

1.28 

2.67 

2.47 

DNR 

10.38 

4.69 

18.11 

25.72 

21 

1.83 

2.76 

2.46 

DNR 

20.35 

7.96 

32.60 

61.39 

22 

1.26 

2.22 

2.01 

DNR 

9.97 

4.60 

17.91 

24.84 

23 

1.67 

3.00 

2.74 

DNR 

19.81 

7.91 

32.66 

67.96 

24 

1.52 

3.12 

2.91 

DNR 

11.71 

5.59 

25.27 

21.57 

25 

1.83 

4.17 

3.96 

DNR 

18.27 

8.37 

33.19 

48.40 

26 

1.08 

4.45 

4.05 

DNR 

11.38 

5.51 

25.05 

19.34 

27 

1.62 

4.42 

4.21 

DNR 

16.37 

7.50 

30.45 

41.98 

28 

4.40 

6.35 

5.37 

DNR 

DNR 

13.91 

53.87 

83.98 

29 

4.87 

7.39 

6.25 

DNR 

DNR 

14.51 

52.55 

117.83 

30 

4.88 

9.08 

7.90 

DNR 

DNR 

16.00 

61.33 

152.21 

31 

5.68 

9.59 

7.58 

DNR 

DNR 

17.05 

61.33 

135.73 

32 

7.42 

15.70 

11.73 

DNR 

DNR 

22.88 

78.63 

553.93 

33 

7.82 

20.20 

15.95 

DNR 

DNR 

25.89 

101.92 

210.14 

34 

8.21 

17.10 

13.76 

DNR 

DNR 

25.42 

92.25 

248.16 

35 

8.81 

19.39 

15.87 

DNR 

DNR 

29.96 

DNR 

DNR 

*DNR — Did  not  run 
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of  the  new  labeling  procedures.  (There  has  indeed  been  some  confusion 
introduced  into  the  literature  by  a number  of  recent  studies  whose  com- 
parisons have  not  been  based  on  fundamental  methodological  differences 
in  labeling  procedures,  but  simply  on  differences  in  pivot  strategies, 
not  clearly  identified  as  such.)  Thus,  the  timesin  Table  I,  while 
extremely  fast,  should  not  be  construed  as  the  best  attainable  with  the 
ARC-II  code.  Preliminary  tests  with  candidate  list  strategies,  not  yet 
refined  to  achieve  the  most  effective  trade-offs  with  the  new  labeling 
procedures,  have  in  fact  resulted  in  times  that  are  roughly  half  of  those 
in  Table  I. 

5 . 3 Memory  Requirements  of  the  Codes 

Table  II  indicates  the  number  of  node  and  arc  length  arrays  required 
in  each  of  the  codes  tested  for  solving  capacitated  problems.  It  should 
be  noted  that  the  program  memory  requirements  of  all  of  the  codes  tested 
were  quite  close  (within  1000  words)  excluding  the  array  requirements. 

Thus  the  important  factor  in  comparing  codes  is  the  number  of  node  and 
arc  length  arrays.  Also,  it  should  be  noted  that  the  primal  and  dual 
simplex  codes  require  one  less  arc  length  array  if  the  problem  is  un- 
capacitated. This  is  not  true  of  the  out-of-kilter  code. 

Since  any  meaningful  network  problem  has  to  have  more  arcs  than  nodes, 
it  is  clear  by  Table  II  that  the  primal  and  dual  simplex  codes  have  a 
distinct  advantage  (in  terms  of  memory  requirements)  over  all  of  the 
other  codes.  Further,  this  advantage  greatly  Increases  as  the  number  of 
arcs  increase  and  if  the  problem  is  uncapacitated.  For  example,  for  a 
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TABLE  II 

CODE  SPECIFICATIONS 


Developer 

Name 

Type 

Number 

of  Arrays 

1. 

Barr,  Glover, 
Klingman 

ARC- I I 

Primal  Simplex  Network 

7N 

+ 3A 

2. 

Barr,  Glover, 
Klingman 

SUPERK 

Out-of-kilter 

4N 

+ 9A 

3. 

Bennington 

BENN 

Non-simplex 

6N 

+ 11A 

4. 

Bray  and  Witzgall 

BSRL 

Out-of-kilter 

6N 

< 

00 

+ 

5. 

Clasen 

SHARE 

Out-of-kilter 

6N 

+ 7 A 

6. 

Glover,  Karney, 
Klingman 

PNET 

Primal  Simplex  Network 

7N 

+ 3A 

7. 

Glover,  Karney, 
Klingman 

DNET 

Dual  Simplex  Network 

9N 

+ 3A 

8. 

Glover,  Karney, 
Klingman 

PNET- I 

Primal  Simplex  Network 

6N 

+ 3A 

9. 

General  Motors 

GM 

Out-of-kilter 

3N 

+ 6A 

N - Node  Length 
A - Arc  Length 
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problem  which  has  10  times  as  many  arcs  as  nodes,  ARC-II,  PNET,  PNET-I, 
or  DNET  require  only  about  one-half  the  memory  space  of  the  best  of  the 
other  codes. 
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